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We use the Schwinger-Dyson equations in the presence of a thermal bath, in order to study chiral 
symmetry breaking in a system of massless Dirac fermions interacting through pseudo quantum 
electrodynamics (PQED3), in (2+1) dimensions. We show that there is a critical temperature T c , 
below which chiral symmetry is broken, and a corresponding mass gap is dynamically generated, 
provided the coupling is above a certain, temperature dependent, critical value a c . The ratio between 
the energy gap and the critical temperature for this model is estimated to be 27r. These results 
are confirmed by analytical and numerical investigations of the Schwinger-Dyson equation for the 
electron. In addition, we calculate the first finite-temperature corrections to the static Coulomb 
interaction. The relevance of this result in the realm of condensed matter systems, like graphene, is 
briefly discussed. 
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I. INTRODUCTION 

Chirality is a physical quantity, carried by fermionic 
particles, which characterizes how these behave under a 
specular reflection. Chiral symmetry is the invariance 
under an operation relating different chiralities. It may 
be either discrete or continuous but, in any case, mass¬ 
lessness is usually a necessary condition for a system to 
be chiral symmetry invariant. Chiral symmetry is usu¬ 
ally broken at a quantum level by radiative corrections, 
which are non-perturbative in the coupling constant. 

A very convenient platform to investigate chiral sym¬ 
metry breaking is the non-perturbative approach pro¬ 
vided by the Schwinger-Dyson equations (SDE). This ap¬ 
proach consists of a complicated system of integral equa¬ 
tions relating exact Green’s functions and vertex func¬ 
tions pj. Nontrivial solutions of the SDE usually im¬ 
ply a dynamical mass generation for the matter field. 
The case of quantum electrodynamics in (3+1) dimen¬ 
sions (QED4) has been studied in Ref. [2j. For quantum 
electrodynamics in (2+1) dimensions (QED3), the chiral- 
symmetry breaking for massless Dirac fermions has been 
extensively studied, both at zero 0 , 3 |t> and finite tem¬ 
perature [MS- 

Three-dimensional theories, such as QED3, have re¬ 
cently attracted much attention because these models are 
relevant for the theoretical description of some effects 
typically observed in condensed-matter systems. Ex¬ 
amples range from high-T c superconductivity 00 to 
graphene [12H16| . Graphene is particularly interesting be¬ 
cause its linear tight-binding dispersion relation coincides 


with the one of massless Dirac particles. Furthermore, 
the system has particle-hole symmetry 0 - 

However, even in two-dimensional materials the real 
photons propagate in (3+1) dimensions. This raises the 
question of how to describe the electromagnetic inter¬ 
actions of charged particles in _{2+l) dimensions. This 
problem was studied in Refs. 0-0. After projecting 
the photons onto (2+1) dimensions, a non-local term 
emerges in the Maxwell Lagrangian, leading to the so- 
called pseudo quantum electrodynamics (PQED) 0 , 0 , 
sometimes also referred to as reduced quantum electro¬ 
dynamics 0. frit 

It has been shown that PQED correctly describes the 
physical 1/r Coulomb interaction between static charges, 
contrarily to QED3 at tree level, which provides a loga¬ 
rithmic static interaction ln(e 2 r). The processes of canon¬ 
ical quantization of theories such as PQED may be found 
in the Refs [24l - l26| . There, it was shown that PQED re¬ 
spects causality. Furthermore, it has been shown that 
PQED respects the Huygens principle in (2+1) dimen¬ 
sions, unlike QED3 0]. More recently, it has been also 
proved that this theory respects unitarity 0. 

In a previous publication, we investigated the occur¬ 
rence of chiral-symmetry breaking in PQED at T=0 0]. 
In this paper, we study the effect of temperature on the 
dynamical chiral-symmetry breaking in PQED coupled 
to massless Dirac fermions, by considering this theory in 
the presence of a thermal bath. At the classical level, the 
theory is chiral invariant due to the absence of a mass 
term for the fermions. We show that there exists a criti¬ 
cal temperature T c , below which quantum effects produce 
dynamical breakdown of the chiral symmetry, with the 
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associated generation of a mass gap. This occurs above 
a critical coupling a c (T), which is temperature depen¬ 
dent and increases with temperature. The temperature, 
therefore, is an inhibitor of dynamical mass generation. 
We estimate the ratio between the mass gap and T c to be 
approximately 2-7T. All these results are then confirmed 
by numerical calculations. Finally, we evaluate the cor¬ 
rections to the static Coulomb potential of PQED, due 
to the finite temperature. For this purpose, we use the 
SDE for the gauge-field propagator, which involves the 
vacuum polarization. We find a logarithmic correction 
in the limit of short-range distance. For the long-range 
limit, the corrected interaction is proportional to the in¬ 
verse of the third power of the distance. 

The outline of this paper is the following: In Sec. II we 
review the Feynman rules for PQED, while in Sec. Ill we 
study the electron self-energy in PQED by using the SDE 
and the imaginary-time formalism or Matsubara frequen¬ 
cies. In Sec. III.A we prove that the sum over Matsubara 
frequencies is convergent for PQED, thereafter we solve 
this sum; in Sec. III.B we use the zero-mode approxi¬ 
mation in order to solve the SDE analytically, and in 
Sec. III.C we use the zero-external-momentum approach 
in the mass function. We obtain the critical temperature, 
the critical coupling constant, and we also estimate the 
ratio between the energy gap and the critical tempera¬ 
ture. Then, we discuss how to include more Matsubara 
frequencies in the calculations. In Sec. IV we calculate 
the corrections to the Coulomb potential in (2+1) di¬ 
mensions due to the thermal bath. We also include two 
appendices, in the first one we compare our analytical 
results to the numerical solution of the full integral equa¬ 
tion for the mass function, while in the second one, we 
discuss the approach about the wave function renormal¬ 
ization. 


II. CHIRAL-SYMMETRY BREAKING IN PQED 

AT T = 0 

The Lagrangian of the PQED in Euclidean space-time 
and for massless Dirac fermions is given by [181] 

where □ is the d’Alembertian operator, e is the electric 
charge of the electron, is the usual field intensity 
tensor of the U(l) gauge field the pseudo electro¬ 
magnetic field, ip is a four-component Dirac field, and 7 ^ 
are rank-4 Dirac matrices. 

The Feynman rules, at zero temperature, imply that 
the gauge-field propagator reads 

A °- W = i <2) 

in the Landau gauge and that the inverse of the fermion 
propagator is 

S of(p) = (3) 


Recently, the authors in Ref. [l6| investigated chiral 
symmetry breaking in PQED at zero temperature using 
the quenched- and the unquenched-rainbow approaches 
to solve the SDE ( see Ref. [l] for a insightful review 
about the SDE). In the quenched-rainbow approach, the 
vacuum polarization is neglected, whereas in the un- 
quenclied one it is considered. Within the first approach, 
chiral-symmetry breaking only occurs if a > a c , whereas 
for the second, a value for the critical number of fermions 
assets given by N c (a) was obtained. In this case, chiral 
symmetry breaking only occurs if Nf < N c (a). Fur¬ 
thermore, it has been shown that the model has intrin¬ 
sic characteristics of both, QED4 (renormalizable theory) 
and QED3. Specially, for QED3 the critical number of 
active fermions is reobtained from PQED in the strong 
coupling limit, i.e, N c (a —> 00 ) = N® Em = ^/f 2 . For 
QED3 this value was first obtained in Ref. Q. 

III. QUENCHED-RAINBOW APPROXIMATION 
AT FINITE TEMPERATURES 

The chiral-symmetry breaking in QED4 at finite tem¬ 
peratures was obtained from non-trivial analytical solu¬ 
tions for the mass function E(p), which only appear be¬ 
low a certain critical temperature T c , estimated using the 
zero mode of the Matsubara frequencies Q. Above the 
critical temperature, the chiral symmetry is restored. It 
has been also shown that the integral equation for the 
mass function in QED4 accounting for temperature ef¬ 
fects is the same as in QED3 by taking T = 0 and using 
a dimensional coupling constant 7 = e 2 T/( 47 r) = aT, 
fixed at the limit T —> 0. In addition, this analytical 
approach for QED4 has been also confirmed by a numer¬ 
ical solution of the full Schwinger-Dyson equation for the 
mass function [7j. However, the same set of approxima¬ 
tions is not applicable to QED3 because there are infrared 
divergences in the electron self-energy if the vacuum po¬ 
larization is not taken into account. 

Due to the similarities between PQED and QED4, here 
we use the same set of approaches proposed in Ref. 
and as expected, no infrared divergence is observed in 
the SDE. The inverse of the full electron propagator then 
reads 

Sf 1 ip) = S 0F ip) - E ip) > ( 4 ) 

where the electron self-energy S(p) is given by 
r + 3 1 , 

S(p) = e 2 y —^S F {k)Y^{k,p)^ v (p-k), (5) 

with r"(fc,p) denoting the full vertex function. In 
the lowest order, we consider T 1 / (/c,p) = 7 ", which is 
called rainbow approach. Furthermore, we use only the 
free gauge-field propagator. This assumption is called 
quenched approach. Both approximations are moti¬ 
vated because of perturbation theory, which establishes 

r"(fc,p) = Y + 0 (e 3 ) and A M „ = A 0/1 „ + 0 (e 2 ) [11-0] . 
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The inverse of the full fermion propagator can be writ¬ 
ten as jl]n:4S 

S^ip) = -P M 7 M A(p) + £ (p), (6) 


(2 n+ 1 )ttT and ui m = (2m + 1)7 tT, with n and m integers 
[29( 1. The integral over k 0 becomes a sum in n by replac¬ 
ing f dkof(ko, k, T ) —> 2irT /(w„, k, T) and the mass 
function is rewritten as 


where A(p ) is the identity-proportional term (this term 
is also called the renormalization function), and E (p) is 
called the mass function. To obtain £(p) we use Eq. © 
in Eq. Q and calculate the trace over the Dirac matrices. 
The corresponding result for zero temperature is fl6| 


£ m (p) = 


x 


y? r d 2 k 

(2^) 2 (2n + l) 
1 


[4 (n — m) 2 n 2 T 2 + (p 


(47raT)E(fc) 

2 n 2 T 2 + k 2 + £(fc) 2 


k) 2 ] 1 / 2 ’ 


( 12 ) 


E (p) = Ana J 


d 3 k £(fc) 1 

(27r) 3 A 2 (k)k 2 + £(/fc) 2 | p - k \’ 


(7) 


with a = e 2 /47r denoting the coupling constant. To ob¬ 
tain A(p ), we multiply Eq. (U) by the Dirac matrices and 
calculate the trace over them. The result for zero tem¬ 
perature is HI 


A{p) = 1 + 


x 


47ra f d 3 k A(k)A(q) 

~p r J (2tt) 3 k 2 A 2 (k) + E 2 (fc) 

( P-q)(k-q ) 


where 


( 8 ) 


q = p — k 


(9) 


and 


A (q) = 


1 

2(g 2 ) 1 / 2 ' 


( 10 ) 


It has been shown that A(p) ss 1 for zero temperature. 
Therefore, if A(p) = 1, E (jp) fully describes the electron 
self-energy S(p) jl6| . In appendix B, we show that this 
result also holds for finite temperatures. Therefore, from 
now on we consider A(T, p) = 1. 

Eq. O is very similar to Eq. (14) in Ref. [3. However, 
there is a different coefficient in the mass function due to 
a constant factor in the gauge-field propagator. Although 
this will change the critical behavior of the theory, one 
verifies that the correct factor in Eq. © should be 1/2 in 
order to obtain the precise Coulomb interaction between 
charges in the plane. Indeed, the critical coupling con¬ 
stant for the onset of chiral-symmetry breaking in Eq. 0 
is a c = 7r/8 ~ 0.40 for T = 0. We rewrite Eq. 0 as two 
integrals, one for the energy component and another for 
the two-dimensional momentum components, 


E (p) = Ana 


dk 0 f d 2 k E (k) 

( 2 n)J ( 2 ir ) 2 k 2 + k 2 + £(fc ) 2 
1 

X [{Po — ko) 2 + (p — k) 2 ] 1 / 2 ’ (11) 


This is the integral equation for the mass function in 
PQED3 at finite temperatures when considering the sum 
over all Matsubara frequencies in the quenched-rainbow 
approach. We have an additional complication with re¬ 
spect to the zero temperature case, which is the depen¬ 
dence of the E m (p) function on the Matsubara frequen¬ 
cies aim- We will assume that the mass function is inde¬ 
pendent of u) m , and therefore consider m = 0. This also 
implies that £(fc) does not participate on the sum over 
n. 


A. The Sum Over Matsubara Frequencies in PQED 

In QED4, the sum over all Matsubara frequencies 
present in the integral equation for the mass function 
£(p, T) is calculated by standard methods [?j|. However, 
for PQED it is not simple to obtain the solution of the 
complete series. First, we need to analyze its conver¬ 
gence. Starting from Eq. m, we see that the sum over 
n can be expressed as a sum U, given by 

+oo -1-00 1 J 

u - Un= 2 Vn 2 + A 2 (2 n + l) 2 + B 2 ' 

n=—oo n =—oo v v y 

Note that the terms of the sum are all positive, 
since A 2 ,B 2 > 0. Furthermore, we can write U = 
2 J2n=i + u 0 ; where = (u n + U- n )/ 2. Therefore, 
it is clear that if u' r t' > is a convergent sum, then U 

will be too. 

Next, we define another independent sum V, which is 

+oo +oo i i 

v = Y,' )n = Y.\(n + A) (n + BY (14) 

where V is convergent and has a known result. Further¬ 
more, we have that v n and u ^ satisfy the inequalities 
v n > > 0 and 

M C+) 

lim —— = 0. (15) 

n—¥oo v n 


Next, we introduce the effects of temperature through 
the Matsubara frequencies in the temporal component. 
For this purpose, we must replace k^ = (ko, k) —► (u n , k) 
and pp = (po,p) —> (w m ,p), where co n and u> m are 
the Matsubara frequencies for fermions, given by w n = 


Therefore, we can use the comparison test of series [3C|, 
which establishes that if V is convergent, then 
is necessarily convergent, which implies that U is also 
convergent, completing our proof. Furthermore, we have 
uo > u n , for any n. 
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From now on, we will not write time components any¬ 
more, hence it is needless to use bold in the momentum 
components. This allows us to simplify the notation, 

E 0 (p,r) = E(p). 

After showing that U is convergent, we follow an ana¬ 
lytical procedure to calculate the sum over the Matsub- 
ara frequencies in Eq. m- This is useful in order to 
identify the temperature-independent term of Eq. m, 
which is called the vacuum term [3lj . Thereby, we rewrite 
Eq. (fl2l) for m = 0 

/ d 2 k r + °° fill 

/ -%(p, £(*)), (16) 

J^ 7T 

where the kernel a y is given by 

“hOO rj~t 

° v = 4n 2 7r 2 T 2 + uj 2 , (2 n + 1) 2 t r 2 T 2 + e 2 ’ 

n= —oo Pi K x K 

with 

= y 2 + (p~ k) 2 , (18) 


and 


e 2 k = k 2 + T 2 {k). (19) 


where 


G = 




2\/q 2 - e 



(26) 


Admitting E(fc) is finite for k —> oo, it follows that 
Eq. m is logarithmically divergent. This is the same 
degree of divergence of the mass-proportional term in 
the electron self-energy, obtained through perturbation 
theory at one-loop expansion at zero temperature. Here, 
however, since we are interested in the non-perturbative 
regime, we shall consider E(fc) a non-trivial function of 
the momentum, which vanishes at large momentum. For 
this purpose, it is more convenient to convert the inte¬ 
gral equation into a differential equation with asymptotic 
conditions. 

Although the previous analytical method is sufficient 
to solve the sum over all Matsubara frequencies, the an¬ 
gular integration obtained from Eq. US and Eq. ((201) 
does not allow us to obtain the differential equation for 
the mass function. Next, we shall explore a different an¬ 
alytical approach to circumvent this problem. 

We return to the integral equation, but we integrate 
out the angular variable instead of summing over the 
Matsubara frequencies. Therefore, Eq. (El reads 


Next, we calculate the sum over the Matsubara fre¬ 
quencies in Eq. o and find 


-i rri f- OO 

S (P) = ^^/ kdkZ(k)I(p,k,i:(k)), (27) 

( 27r ) Jo 


(jy = C B [2n B (Up tk ) + l]-C F [l-2n F (e k )], (20) where I is given by 


where n B (oJ p ,k ) = [exp (uj Ptk /T) - 1] 1 and n F (e k ) = 
[exp(efc/T) + l] -1 are the Bose and Fermi distribution 
functions, respectively. The functions C B and C F are 


C b (T) = 


{sp,kd Pt k + 7 t 2 T 2 ) 


2 Wp,fc(d 2 jfc + 7r 2 T 2 )(s 2 fc + tt 2 T 2 ) ’ 


and 


C f (T) = 


{s p ,kdp,k - 7t 2 T 2 ) 

2e k (dl k +n^)(slk + n 2 T^y 


( 21 ) 


( 22 ) 


where d Pjfc = e k - w Pjfc and s Pifc = e k + ui Ptk . 

We obtain the vacuum term by using T = 0 in Eq. ©■ 
Thereby, the kernel at zero temperature is 


a v (T = 0) = C B {0)-C F (0) 


1 

2 kJp,k€k(^p,k T C/ c ) 


(23) 


I 

/(fc,p, E(fc)) = (2n + 1) 2 7 t 2 T 2 + k 2 + E(fc) 2 

x K ( x ”) (28) 
\J(jp— k ) 2 + 4n 2 7r 2 T 2 ’ 

with K[x n ) denoting the elliptic function. For x n < 1, 
we have K(x n ) = 7t/2 + irx n /8 + 0(x^), with x n given 

by 


4 pk 

(p — k ) 2 + 4n 2 7r 2 T 2 


(29) 


Next, we propose an approximation to the sum U, in 
order to obtain an analytical solution to £(p). The zero¬ 
mode approximation, i.e, U « uq has been used for QED4 
in Ref. [7] and is also valid for the PQED. 


Integrating out y in Eq. El yields 


2efc-\/q 2 - e; 


1-tan 


-l 


Cfc 




where q = p — k. Therefore, Eq. El becomes 


B. Zero-Mode Approximation 


(24) 


For low temperatures, the sum over n in Eq. (1281) is ap¬ 
proximately given by the term uq. Therefore, by taking 
n = 0 in Eq. m, we find 


S(p)| T =o 


= 47ra 


/ 


d 2 k E(jfc) 

(2tt) 2 k 2 + E 2 (fc) 


G(q,S(fc)), 


(25) 




167T7 r°° k K(x 0 )T,(k) 

(2tt) 2 J 0 |p - k\ tt 2 T 2 + k 2 + E(fc) 2 ’ 


(30) 
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where 7 = aT is a dimensional coupling constant. This 
condition is necessary to study the dynamical mass gen¬ 
eration in the case of a nonzero critical temperature T c . 
However, one may always return to the coupling constant 
and consider it in the context of a finite critical coupling 
constant a c (T) at finite temperatures. The symmetry 
is fully restored if T c = 0 or a c (T) —> 00 . Note that 
both interpretations are equivalent, in the sense that at 
very low temperatures, interaction effects become more 
relevant. 

Next, we divide Eq. (1301) into two parts, the infrared 
part k < p and the ultraviolet part k > p. Furthermore, 
we include a ultraviolet “cutoff” A (bear in mind that 
we may return to the continuum limit A = 00 at any 
time in the calculations). However, for comparison with 
lattice field theory, it is interesting to keep the cutoff 
finite. Using these approaches in the elliptic function, we 
obtain 

^6{p-k) + ^0{k-p). (31) 

Therefore, Eq. becomes 


+2l l ‘' W+^+SM 


2 ' 


(32) 


We can neglect the nonlinear term and approximate 
E 2 (fc) + k 2 + 7 r 2 T 2 ss k 2 + tt 2 T 2 . As argued earlier in 
the literature, this approximation does not change sig¬ 
nificantly the critical behavior of the theory [S]. Further¬ 
more, T is dominant in the region where k —t 0 and is 
negligible in the region where k —> A. This shows that 
the temperature behaves like a natural infrared cutoff for 
Eq. (1551) . This infrared cutoff is identified as the critical 
temperature for chiral symmetry breaking [?}. Thereby, 
by taking the first derivative of the above expression, we 
have 


<ffi(p) 2 7 fP E(fc) 

-nr = -?U ik — 


(33) 


which can be converted to a second-order differential 
equation given by 


d_ / 2 dS{p) \ 
dp V dp ) 


2 7 


£(p) = 0 , 


(34) 


with infrared and ultraviolet boundary conditions given, 
respectively, by 


lim p 2 

p—yrrT c 


d£(p) 

dp 


= 0 


(35) 


and 


lim E (p) = 0. (36) 

p—> A 


The solution of the differential Eq. (1341) is 


£ (p) = c x 





(37) 


where J\(x) and Yi(x) are the Bessel functions of first 
and second kind, respectively. C\ and Ci are arbitrary 
constants with dimension [C 1 ] = [C 2 ] = 1, both in units 
of energy. From Eq. (l36l) . it is possible to show that C 2 = 
0, because xYi(x) is not zero if a; —>• 0, with x = -^Ay/A. 
On the other hand, Eq. (1551) provides the identity 


2Jr(0+^o(0-^2(0=0, (38) 

where £ = ^/ 87 / 7 tT , c . Eq. (1551) has a set of nontrivial 
solutions {£„} = {£ 0 , £ 1 , •■•}, where < £ n+ i for any 
n. The minimal value £0 is chosen to define the critical 
temperature T c , 


T r = 


87 

<o' 


(39) 


Using Eq. (1551) . one can show that £0 ~ 2.40, thus 
allowing us to determine the critical temperature only 
in terms of the dimensional coupling constant, T c = 
8aT/nt;Q = 0.44 aT, i.e, for T < T c , then a > a c (T ss 
T c ) = tt£o /8 ~ 2.26. This result shows that even for 
low temperatures the critical coupling constant is much 
larger than in the case of zero temperature aJT = 0) = 
77/8^0.40 0. 

For an arbitrary value of 7 , it can be verified by numeri¬ 
cal tests that the full integral equation, given by Eq. m , 
obeys E (p,T > T c ) —> 0. Therefore, a nontrivial solution 
for E (p) requires a nontrivial value for T c . 

The approach of this section has provided an analyti¬ 
cal solution for the mass function E (p,T), as well as an 
estimative for T c . Unfortunately, from the analytical so¬ 
lution, it is not clear that E (p,T = T c ) — > 0. In order 
to confirm the phase transition at T c and a c (T), in the 
next section, we will recalculate the mass function, but 
at zero external momentum. 


C. Zero-External-Momentum Approximation 


In this section, we will take zero external momentum 
in Eq. (1501) . i.e, E {p,T) m E(0,T) = m(T). At zero 
momentum, the mass function is the pole of the fermion 
propagator, therefore 2 m(T) is the energy gap between 
positive and negative energies. It has been argued that 
the mass function reaches its maximum value exactly at 
p -A 0 @,0. Therefore, if m(T) / 0 we immediately 
conclude that chiral-symmetry breaking will occur. By 
using this approach into Eq. we have 


1 = 27 


dk 


0 


tt 2 T 2 + k 2 + m(T) 2 ' 


(40) 
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Eq. ® can be transformed into a transcendental 
equation for m(T), given by 


Jm(T ) 2 + 7r 2 T 2 , 

——-= tan" 1 

2y 


A 


sJm(T) 2 + tt 2 T 2 


■ (41) 


In order to obtain an analytical solution for m(T), 
we expand the right-hand side (rhs) of Eq. (HT1) for 
A > m(T) 2 + 7 r 2 T 2 . This is a reasonable assumption, 
since we are in the low-temperature regime and the gen¬ 
erated mass is expected to be much smaller than the 
ultraviolet cutoff. We perform a Taylor expansion until 
second order, tan" 1 (a;) ~ 7r/2 — 1/a: + 0{x ~ 2 ) if x > 1, 
and we find two possible solutions 


m(T) » ±ny/(T2-T*), (42) 


analytical and numerical approaches by using this regu¬ 
lator in the present work. 

We now connect T c with the mass function at zero 
temperature, using Eq. (1421) . For T = 0, we obtain T c « 
m(0)/7T. Thus, it is possible to calculate the ratio R 
between the energy gap 2?n(0) and T c , given by R = 2i r. 
For QED3, the authors in Ref. 0 calculated this ratio 
and obtained Rqed 3 = 9.36 for one copy of the Dirac 
field. 

Eq. (l42l) shows that the chiral-symmetry breaking in 
PQED at finite temperatures is a second-order phase 
transition, with critical exponent equal to 1/2. 

The zero external-momentum approximation is an 
ideal approach to include the sum over Matsubara fre¬ 
quencies in Eq. m■ Indeed, in this limit we have 


where T' 2 = 7 2 A 2 /( 2 y + A) 2 . From now on, let us admit 
that m(T ) > 0. Next, we return 7 — > aT, to obtain the 
critical coupling constant 

»‘< T > = (a< 43 » 

Eq. (l43l) allows us to estimate an interval for a c (T). 
Let us choose T = 0.3A as an upper bound, because 
T/A <C 1. In this case, a c (T = 0) < a c < a c {T = 0.3A), 
i.e, 1.0 < a c < 2.5. For the particular value T = 0.5A 
the chiral symmetry is completely restored. However, it 
is clear that the high-temperature limit T « A —> 00 is 
not allowed within this approach. 

Having in mind condensed-matter applications, we 
consider the fact that the cutoff A is a physical param¬ 
eter determined by the lattice spacing, i.e, A oc 1/a^, 
where usually a l ~ 10" lo m. For graphene, the cutoff is 
A = Hvf/cil ~ 3.0 eV [Hj]. Conversely, from a quantum 
field theory perspective, one should take A = 00 , which 
implies a c = 1 . 0 . 

In the case of quantum field theory applications to par¬ 
ticle physics, however, dimensional regularization would 
be more appropriate. In this case, we should return 
to Eq. (12) and perform the dimensional regulariza¬ 
tion, generalizing the two-dimensional integral into a D- 
dimensional integral. Then, after some calculations, we 
conclude that the critical coupling in the zero-external- 
momentum approach is equal to one. In fact, this result 
may also be obtained from Eq. flUl) with A = 00 . There¬ 
fore, a c (T) = 1.0 is a result that does not depend on 
the regulator. Nevertheless, the full understanding of 
the pattern of chiral-symmetry breaking in PQED, using 
the dimensional regularization, requires a more complete 
study of the integral equation in Eq. m, especially ap¬ 
plying better analytical and numerical approaches that 
go beyond the scope of the present work. 

A similar problem has been discussed in the frame¬ 
work of quenched quantum electrodynamics in (3+1)- 
dimensions (QED4) §54], 35|. For QED4, it has been 
shown that the results in these two different regulators 
are in agreement. We shall not discuss more elaborated 


+OO 


1 = 2y / dkk 


1 


n =—00 


(2 n + 1) 2 t t 2 T 2 + k 2 + m(T) 2 
1 


y/k 2 + 4n 2 7r 2 T 2 


.(44) 


We can find numerical solutions to m(T) by using 
Eq. (0U), thereby verifying that the zero-mode approx¬ 
imation is in good agreement with the SDE for all Mat¬ 
subara frequencies. This concludes our investigation of 
chiral-symmetry breaking within the quenched-rainbow 
approach. The general important points to be retained 
are: (i) The quenched-rainbow approach provides reason¬ 
able physical solutions for the mass function in PQED, 
similarly to QED4 and quite differently from QED3; (ii) 
Nontrivial solutions have been found for T < T c \ (iii) T c 
is estimated from the energy gap at zero temperature; 
and (iv) The zero-mode approximation is found to be in 
agreement with the sum over all Matsubara frequencies. 

We have assumed that H(p,T > T c ) -+ 0 based on 
the fact that systems at high temperatures naturally are 
less prone to quantum effects. Then, no chiral symmetry 
breaking should be observed. However, for the sake of 
completeness, we will briefly discuss how to obtain ex¬ 
plicitly this trivial solution. From the integral Eq. CP, 
we must use p 2 + T 2 + E 2 (p) ss T 2 , a realistic assumption 
for high temperatures. Then, by using the expansion in 
Eq. (I3T1) . we find the differential equation 


d_ / 2 dS(p) \ 
dp V dp ) 


= 0 , 


(45) 


the general solution of which is 

£(p) = — + C 4 . (46) 

P 


Note that C\ = 0, since E(p -+ 00 ) = 0, and C 3 = 0 
because E(p —> 0) must be finite. In this limit, T is no 
longer an infrared cutoff, and T,(p,T) = 0. 
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IV. STATIC INTERACTIONS BETWEEN 
MASSLESS DIRAC PARTICLES AT T^O IN THE 
PQED APPROACH 


where y = pr. The first term in the rhs of Eq. (l55j) is the 
Coulomb interaction, while the second one is a correction 
due to the thermal bath. The integral over y yields [32j 


In this section, we recalculate the static interaction be¬ 
tween charges in the plane, in order to clarify the phys¬ 
ical interpretation of the electromagnetic interaction in 
the plane at finite temperatures. For this purpose, we 
must include the corrections of the vacuum polarization 
to the gauge-field propagator. 

From Eq. 0 , we obtain the static gauge-field propa¬ 
gator, given by 

A^(p 0 = 0,p) = ( 47 ) 

It has been shown [l 8 j that this propagator reproduces 
the Coulomb interactions for static charges in the plane. 
The static vacuum polarization at finite temperatures has 
been calculated in the Ref. 51 , and for massless Dirac 
particles it is given by 

2 a T 

iWpo = o,p,r) =-V><W(p,r), (48) 

7r 

where 



dy 


My) 

2 y + rM 2 (T) 


7r 

4 





,(56) 


where 


ro(T) 


2 

M 2 (T) 


7r 

a T In 2 ’ 


(57) 


H 0 (r/r 0 ) is the Struve function, and Y 0 (r/ro) is the 
Bessel function of the second kind. Therefore, the static 
interaction between charged particles in PQED in the 
presence of the thermal bath is 


V(r, r 0 ) 



. (58) 

We now rewrite the potential V(r,ro) as a function of 
the dimensionless variable l = r/ro(T ), 


font) = ^i{ l -\ l l H o (0 - >0(1)]} ■ (59) 


f(.P,T)= [ dxhi 
Jo 

Let us now to consider the high-temperature limit p 
in Eq. (l49l) : then /(p —> 0,T) = In2. 

The corrected propagator reads 

A^(p,T) = A"J(p) - n M „(p,T), 



(49) 

-A 0 


(50) 


For l -C 1, we have 


r 0 V(l) 



1 + Zln 




(60) 


where In 7 = — y e and y e is the Euler’s constant. For 
l 1, we use the asymptotic representation of the Struve 
function (32[ 


leading to 


where 


A /lt /(p, T) = 




2y / p 2 + M 2 (T) ’ 


M 2 (T) = 


2 aT In 2 


(51) 


(52) 


(61) 

(62) 


and consider p = 2 , to find 


rom = ^\^+°(i- 2 )]. 


The static interaction V(r) is obtained from a Fourier 
transformation of Eq. m, 


V (r) = e 2 


d 2 p e~ ipr 
(2tt) 2 2Vp 2 + M 2 (T)' 


(53) 


The angular integral in Eq. (l53ll is performed, yielding 


V(r) 


2 f °° dp pMpt) 

J Q 2 tt 2 p + M 2 (Ty 


(54) 


where Jo(pr) is the Bessel function and p = y/p 2 . By 
defining z = 2p + M 2 (T) and changing the integra¬ 
tion variable into z, we obtain, after some algebra, that 
Eq. (15H) may be rewritten as 


M 2 (T) 


V(r) = - -e 2 —- 

47rr 47 t 


dy 


Mv) 


2y + rM 2 (T) ’ 


(55) 


The last result indicates that the corrected potential, 
for r ro(T), goes to zero much faster than the usual 
Coulomb potential. In Fig. 1, we plot Eq. (l58l) and the 
usual Coulomb potential for comparison. 

Next, we estimate the physical value of r 0 (T). Let us 
assume a massless particle with energy dispersion E = 
vp = hv/r, where v = fc is the particle velocity, 0 < / < 
1, and c is the light velocity. For graphene / = 1/300 
and for a relativistic particle / = 1, for instance. Taking 
this into account, we rewrite Eq. (1571) in physical units, 


ro(T) = 


nHfc 

af IcbT In 2 ’ 


(63) 


where we reintroduced h and kg, the Planck and Boltz¬ 
mann’s constant, respectively and af = aQEm/f = 
1/(137/) is the coupling constant. Thus, ro(T) ~ 


























FIG. 1: (Color online) We artificially consider a = 1/137 
and ro = 1.0 for both lines. The thick line is the corrected 
potential given in Eq. (1581) . while the dashed line is the usual 
Coulomb potential for T = 0. Here r and r o are given in units 
of T" 1 . 

10~ 5 /T m, with T in Kelvin and / = 1/300. At room 
temperature T = 300 K, we find ro ss 10 nm, which is 
much larger than the atomic scale. This value is smaller 
when / is less than 1/300. 

There is a parallel between the static interaction in 
PQED at finite temperatures and QED3 at zero temper¬ 
ature. Indeed, starting from the screened static gauge- 
field propagator in QED3, obtained from the respective 
SDE, it has been shown that at large r the static electron- 
electron interaction is given by V(r) oc ln(e 2 r) + h(r), 
where h(r ) oc 1/r + 0(l/r 2 ) |3(J. Therefore, the results 
are similar to PQED. However, for QED3 the Coulomb 
potential is the correction to the confining potential and 
the term which breaks scale invariance is the dimensional 
electric charge, instead of the temperature as in PQED. 
For finite temperatures , the sta tic potential of QED3 be¬ 
haves like V(r) oc 1/y/M ( T)r exp(—M ( T)r) [lj|, which 
is quite different from the corresponding result for PQED 
in Eq. (16U1) . 


V. DISCUSSION 

In this paper, we show that the inclusion of finite- 
temperature effects into the SDE leads to a critical pa¬ 
rameter T c for the onset of dynamical mass generation. 
Moreover, at high temperatures the chiral symmetry is 
restored. The finite temperature also increases the value 
of the critical coupling constant; for low temperatures, 
we estimate that 1.0 < a c (T) < 2.5, which are large val¬ 
ues compared to the equivalent zero-temperature result 
a c (T = 0) = 7t/8 ss 0.40. Our analytical findings were 
verified by performing numerical tests. We estimate the 
ratio between the energy gap and the critical tempera¬ 
ture for PQED to be R = 2i r. Furthermore, we show 
that the Coulomb potential between static charges in the 
plane is corrected by a logarithmic potential proportional 


to M 2 (T) oc aT, where a = e 2 /( 47 r) for the short-range 
limit, while a third-power potential emerges in the long- 
range limit. A deeper investigation of the ground state 
generated by this interaction may bring interesting phys¬ 
ical results. 

For graphene, the investigation of chiral symmetry 
breaking in ( 2 + 1 ) dimensional theories is an important 
topic, since this effect has been related to a dynamical 
gap generation, which can be relevant to describe elec¬ 
tronic transport in this system. It has been argued that 
this gap energy leads to important technological appli¬ 
cations for graphene, such as a graphene transistors |37| . 
PQED at T = 0 admits gap generation, but only above a 
critical coupling constant. Although our results include 
finite temperatures, and thus are more realistic because 
any measurement of the gap is done at T ^ 0, we are 
considering vf = c. The inclusion of still more physical 
parameters to study chiral symmetry breaking in PQED 
may lead to results closer to the experimental findings 
in graphene. It may as well bring a deeper understand¬ 
ing of how to control this energy gap. These parameters, 
among others, should be: the Fermi velocity vf , finite 
temperature T, magnetic field B, and chemical potential 
/ 1 . A recent experimental measurement at low magnetic 
fields has determined an upper limit of 0.1 meV for a 
possible gap in suspended monolayer graphene [Hj. 

According to the Coleman theorem, there is no Gold- 
stone boson in ( 1 + 1 ) dimensions, i.e, no spontaneous 
(continuos) symmetry breaking may occur [39l - l41| . Be¬ 
cause of this theorem, one may conclude that there is no 
ferromagnetism or antiferromagnetism in the Heisenberg 
model in one or two dimensions [icj. For (2+1) dimen¬ 
sional theories at finite temperatures, the theorem also 
applies due to the loop integrals, which behave effectively 
in (1+1) dimensions PI- This is observed for the case 
of PQED in Eq. (1121) . Therefore, there is no continuous 
chiral-symmetry breaking from the very beginning and 
the dynamically generated mass instead breaks the dis¬ 
crete chiral symmetry. For the case of the continuous 
chiral symmetry, the dynamical mass generation could 
be related to the Kosterlitz-Thouless mechanism [HI, as 
has been discussed for the Gross-Neveu model in (2+1) 
dimensions at finite temperatures [Hj. Because of the 
Kosterlitz-Thouless mechanism, no continuous symmetry 
is broken, even in the massive phase, thus avoiding any 
possible disagreement with the Coleman theorem [43[ . 

The non-local approach of PQED has been applied to 
the study of superconductivity in ( 2 + 1 ) dimensions in 
the view point of the Kosterlitz-Thouless mechanism 0] ■ 
In this case, because of the coupling constant between 
the matter and the bosonic excitation (lattice effect), it 
has been shown that the gauge-field propagator has finite 
mass, thus finite penetration depth, which may lead to 
an effective description of the Meissner effect in super¬ 
conductors. Furthermore, we have obtained the ratio R 
between the energy gap and the critical temperature for 
the matter field, which is almost twice the value obtained 
from the BCS theory for superconductivity. We shall ex- 
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plore elsewhere the possible relation between dynamical 
mass generation, Kosterlitz-Thouless mechanism, and su¬ 
perconductivity [20]. 
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VII. APPENDIX A: NUMERICAL RESULTS 
FOR THE MASS FUNCTION 


In this appendix, we perform some numerical tests to 
verify the validity of the analytical approaches adopted 
in this paper. For simplicity, we will begin by studying 
Eq. m and its analytical solution, Eq. (1421) . in which for 
A large, we obtained T c = 7. In Fig. O we plot the nu¬ 
merical results for A = 10 and 7 = 1 / (87t) with dots, and 
the analytical result as a solid line. Numerical solutions 
are found by looking for roots of the transcendental equa¬ 
tion (l4Tjl . A very good agreement is observed. Note that 
both solutions vanish almost at the same point T c = 0.04 
in units of A/10 (u.A). 



FIG. 2: The solid line is the analytical solution given by 
Eq. (1421) and the dots are the numerical results. We use 
A = 10 and 7 = l/(87r) u.A. The critical temperature is 
T c = 7 = 0.04 u.A. 


Next, we investigated the influence of all Matsubara 
frequencies for the mass function. After calculating the 


integral over the momentum k in Eq. (1441) . we find 

+ OO 


F{m) = 1-27 


1 


tan 


tan 


1 /m 2 + (1 + 4 n) 7 r 2 T 2 

-! C VA 2 + 4 n 2 7 T 2 T 2 ^ 

\ sj m 2 + (1 + 4 n) 7 r 2 T 2 ) 

_i ( 2mrT \ 

1 \Jnri 2 + (1 + 4 n) 7 r 2 T 2 ) 


(64) 



FIG. 3: The dashed line is the numerical solution of the non¬ 
linear integral equation given by Eq. J30J. The solid line is the 
analytical solution Eq. 071) for A = 0.15 and B = 0. For these 
two lines we used A = 10 and T = 0.01 u.A. The thick line 
is the numerical solution for T = 0.03 u.A. For this approach 
T c = 0.447 = 0.018 u.A for 7 = 1/(8tt) u.A. 

The roots of F{m{T)) = 0 yield the mass function 
m{T) for a given value of T. We estimate the critical tem¬ 
perature by imposing that for T = T C the mass vanishes, 
m(T = T c ) = 0. Thus, we obtain a critical temperature 
T c = 0.10 u.A. This result shows that from the very sim¬ 
ple zero-mode approach is already enough to provide an 
estimate of T c in units of A. However, a more accurate 
result is expected when one includes all the Matsubara 
frequencies. To investigate this more complicated case, 
we include an effective limit to the Matsubara frequen¬ 
cies, n € [—Iot+Io] with lo = 10. It is possible to show 
that for any number large than Iq, the function F(m(T)) 
remains nearly the same. Indeed, as we demonstrated, 
the sum over n is convergent, therefore we expect that 
u n = 0 if n —> 00 . 

For non-zero external momentum, we can compare the 
analytical solution in Eq. (1371) to the numerical solutions 
of the nonlinear integral Eq. m ■ For this purpose, we 
have to choose a value for the arbitrary constant A in 
the analytical solution. For temperatures larger than 
T c = 0.447 = 0.018 u.A for 7 = 1/(87 t) u.A, the mass 
function given by the nonlinear integral equation van¬ 
ishes, as expected from the analytical solution. There¬ 
fore, Fig. 0 shows that the analytical solution found ear¬ 
lier is in very good agreement with the nonlinear integral 
equation. 
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VIII. APPENDIX B: THE RENORMALIZATION 
FUNCTION A(p) 


In this appendix, we show that A(T, p) = 1 is also 
a suitable approach for finite temperatures. This ap¬ 
proach has been discussed for QED3 at zero tempera¬ 
ture in Ref. Q, while for PQED it has been verified in 
Ref. [l6j]. After introducing the Matsubara frequencies 
in Eq. © as in Sec. Ill, we obtain 


A(p) = 1 


AiraT 



d 2 k A(k)A(q) 

(2tt) 2 fc 2 A 2 (fc) + E 2 (fc) 


x 


(p.q)(k.q) 


(65) 


where p = (nT, p), k = ((2 n + 1)7tT, k), and q = 
(—2mrT,q). We assume m = 0 because A(p) = A(T, p) 
should be independent on Matsubara frequencies. For 
simplicity, let us consider y/p 1 2 = p and p.k = pk cos 9. 
The angular integral may be solved analytically, and af¬ 
ter some simplifications we obtain 


Jo 9 

8nn 2 T 2 H n (p, k)Y n (p , - k ) 


( pk ) 3 4 5 6 / 2 


E(-4Y n (p,k)), (66) 


where E(x) = tt/2—(tt/8)x+... for \x\ < 1 is the complete 
elliptic integral, 


and 


Yn{p,k) 


pk 

(p — k ) 2 + ( 2nirT ) 2 ’ 


(67) 


H n (p, k ) = 2nY n (p, k) 1/2 (p 2 - n 2 T 2 ) 

+2 pYY 2 (p, k)(p - k) - pkY~ 1/2 (p, k). (68) 

Next, we obtain the integral equation for A(p) with a 
kernel only dependent on the internal-momentum inte¬ 
gral, 


A(p) = 1 


aT I" kdkA(k)I^(p, k) 

J fc 2 A 2 (fc) + E 2 (fc)' 


As expected, A(T, p) = 1 + 0(a), which indicates the 
accuracy of the approximation A(p) « 1 [E]. Beyond this 
argument, we need to look for a numerical solution of 
Eq. (IMll for Y(k,T) = 0 (the symmetric phase), as in 
Fig. 3. Note that, since Y,(k,T) oc a, we should not con¬ 
sider E(fc,T) ^ 0 in order to solve the integral equation 
for A(T , p) because these new terms are of order a 2 . This 
symmetric phase for A(T = 0, p) has been also discussed 
in the context of QED3 [E|. 

In order to obtain the results in Fig. 1 and Fig. 3, we 
converted the momentum-dependent kernel in Eq. (13011 
1.2 F 

1 . - 

0.8 - 

P 

q10.6 - 

< 

0.4 - 
0.2 - 

o L_i_i_i_i_i— 

2 4 6 8 10 

P 

FIG. 4: (Color online) Numerical solution of Eq. (l69l) for 
a — 1/(87r), T = 1.0 u.A, A = 10, n £ [-10,-1-10], which 
shows that A(p) « 1 is indeed a realistic approximation. 


and Eq. (16911 into a system of nonlinear algebraic equa¬ 
tions by using the repeated trapezoidal quadrature rule, 
as has been successfully done for PQED at T = 0 
[r|. The fundamental step is to replace the contin¬ 
uum variables p and k into a set of discrete variables 
Xi and i/i, with 1CD 3 < Xi,yi < A = 10. In this 
interval, we obtained a set of M = 300 solutions for 
the mass function Ei(xj), which implies in a separation 
h = (10 — 10” 3 )/(300 — 1) ss 0.033 between the discrete 
variables. For T / 0, the temperature is only a new pa¬ 
rameter inside the kernel of the nonlinear integral equa¬ 
tions. More details about the numerical procedures may 
be found in Ref. [l6j for PQED at T = 0, in Ref. @ for 
QED3 at T = 0, and in Ref. 0 for QED4 at T / 0. 
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